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1. Introduction 

The two items in the title are rooted in disjoint areas. The Airy process naturally 
belongs to random matrix theory, while the one-dimensional Kardar-Parisi-Zhang (KPZ) 
equation is the stochastic evolution for droplet growth in a thin film. The purpose of 
the introduction is to put forward the two players, starting with stochastic growth. 

The KPZ equation [1] describes the evolution of the height function h{x, t) on the 
line and reads 

dth{x,t) = \\{d^h{x,i)f + udlh{x,t) + /Dr/(x,t) . (1.1) 

Here the first term is the nonlinear growth velocity, to the lowest relevant order in 
the slope, and the Laplacian represents surface tension. The randomness of local 
attachments at the interface is modeled by Gaussian white noise, rj, with mean zero 
and covariance 

{r]{x, t)7]{x', t')) = 6{x - x')5{t - t') . (1.2) 

\/15 is the noise strength. 

As an evolution equation, (11. ip has to be supplemented with initial conditions. We 
are interested in droplet growth which is triggered by the sharp wedge initial condition 

/i(x,0) = -|x|/e (1.3) 

in the limit e — )■ 0. The solution to the KPZ equation is then of the form [3 [3l H] 

r'Kx, t) = - ^ + log/x + e(x, t) (1.4) 
with parameters 

Note that the random amplitude C(x,t) is stationary in x for fixed t. The KPZ theory 
predicts $,{x,t) to be of order t^^^ and to have correlations in x over a range of order 

^2/3. 

To proceed to the second player, the Airy process is most naturally introduced 
through Dyson's Brownian motion. Let us consider the dynamics of an x complex 
hermitian matrix, A{t), as governed by the Langevin equation 

j^A{t) = -^A{t) + W{t). (1.6) 

Here W{t) is unitarily invariant N x N matrix- valued Gaussian white noise, i.e. W{t) 
is a.n N X N matrix, W{t) = W{t)*, such that in some fixed basis 

{W,,{t)W,:y{t')) = 5,y5,,5{t - t') . (1.7) 

In particular, for any unitary matrix U it holds 

UW{t)U* = W{t) (1.8) 
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in distribution. The Langevin equation f ll.6p has a stationary solution, which is given 
by the Gaussian unitary ensemble (GUE) 

Z-^e-TT^^^'dA. (1.9) 

If one picks (11. 9p as initial distribution, then A{t) is a stationary matrix- valued stochastic 
process. Some given trajectory in matrix space induces the motion of the eigenvalues of 
A{t). They can be ordered as 

Ai(t) < ... < Ajv(t) , (1.10) 

since with probability one there is no crossing. We are interested in the largest eigenvalue 
AAr(t). It is roughly located at the band edge 2A^. Rescaling relative to the edge, one 
has 

lim N-^''\Xn{N^'^w) - 2N) = A(w) , (1.11) 

A^— !>oo 

which we regard as the definition of the Airy process A{w). Since XN{t) is stationary, 
so is the limit A{w). 

(11.111) is a somewhat indirect prescription. A more direct characterization comes 
from the fact that the set of trajectories {Xj{t),t G M, j = 1,...,A^} form a 
determinantal process. Thus the finite dimensional distributions of A{w) can be written 
as a Fredholm determinant. At fixed w the relevant operator is the Airy operator K 
with the integral kernel 

POO 

K{x,y)= / duAi{u + x)Ai{u + y) , (1.12) 
Jo 

called Airy kernel, where Ai is the standard Airy function. The transfer operator along 
the ti^-axis is governed by the Airy Hamiltonian 



H=-^^u. (1.13) 

Note that K is the projection onto all negative energy eigenstates of H, thus [H, K] = 
and = K. In particular one can write K{x,y) = {Airc\K\Aiy), where Aix{u) = 
Ai{u + x) is the eigenfunction of H with eigenvalue —x. Now let us pick times 
wi < . . . < Wn and real numbers ai, . . . , a„. Then 

FToh{A{wj) < aj,j = 1, . . . ,n) = det(l - K) . (1.14) 

K is a n X n matrix with operator entries, 

U\K\y)-S l{.>M(a:|e("^)^ir|y)l|,>„^,} for z>j, 

^ ' ^^'^^ \ -l{.>a.}(x|e("^)^(l-ir)|y)l|,>„^,} for z<j. ^'-'"^ 

Kij is trace class and the Fredholm determinant (11.141) is well defined. The marginal 
distribution at one single point equals det(l — PaKPa) with Pa the projection operator 
onto the interval [a, oo), which one recognizes as the GUE Tracy-Widom distribution 
function. 
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One knows that the Airy process is almost surely continuous and looks locally like a 
Brownian motion. In particular for the covariance, {{A{w) — ^(0))^) = \w\ for \w\ — 
and 

(^(0MH)-M(0))(^H)^ A (1.16) 

for large \w\. We refer to [5] for more detailed information. For polynuclear growth [6] 
and the single step model [7] it has been established that in the long time limit, on the 
scale t^/^ for the height and t^^^ in the transverse direction, the height statistics is well 
approximated by the Airy process. On this basis one expects the same property to hold 
for the KPZ equation. More specifically, for the fluctuation term in f ll.4p it should hold 

lim r^/^^(2/it2/3w, 2t) = A{w) , (1.17) 

where the nonuniversal coefficient fi, see (11.51) . results from the KPZ scaling theory [Sj. 

For a single point, n = 1, the property (11.17P has recently been estabhshed. The 
results in [21 El IH E] rely on the approximation of the KPZ equation through the weakly 
asymmetric simple exclusion process. Alternatively one can use the replica method 
[TUl [TTl \T2[ [T3] . The m-th replica is of order exp[tm^/3] and one is forced to make an 
analytic continuation in m into the complex plane C. At this point one adopts the 
choice 

exp[tmV3] = [ dnAi(M)e""*'^' , (1.18) 



which is supported by the fact that thereby the one-point generating function agrees 
with the one in [21 13] . 

For multi-point distributions, currently only the replica method seems to be 
available. We regard it as compelling that the analytic continuation should be chosen 
as in (I1.18p . For two points, Wi, W2, such a program is carried out in [14J. Under a 
decoupling approximation it is established that t~^/^^(2/it^/^Wj, 2t), j = 1,2, jointly 
converge to A{wi) , A{w2) as t — )■ oo. In fact, one first computes the generating function 

(exp(-e-^^+«("^'*) - e-^2+«("2'*))) (1.19) 

in terms of a Fredholm determinant, which is sufficiently explicit to easily read off the 
long time limit. 

The goal of this contribution is to extend the results from [T3] to an arbitrary but 
finite number of points Xi, . . . , x„. In particular we will study the generating function 

n 

exp (- J]e-^^+«("^'*))) . (1.20) 
i=i 

In the properly scaled limit t — )■ oo one obtains a Fredholm determinant of an operator 
acting in L^(M), which looks very different from the standard definition given above. In 



Appendix A the equality of both expressions will be established. 



In the next section we recall the replica method based on the link of the KPZ 
equation to the continuum directed polymer. We also state our main result. In Section [3] 
we carry out the replica summation, which is used in Section IHto compute the generating 
function ( OOD . 



The One- dimensional KPZ Equation and the Airy Process 



5 



(Z,,,,^(x,t)) = ^== exp i-ii) . (2.5) 



2. Replica approach and main result 

Through the Cole-Hopf transformation, 

=e^^("'*) , (2.1) 
the KPZ equation becomes 

dtZ{x, t) = v dlZ{x, t) + {XVD/2iy) r]{x, t)Z{x, t) (2.2) 
and the sharp wedge initial condition translates to 

Z{x,0) = S{x). (2.3) 
The solution of fl2.2p . (12. 3p can be written through the Feynman-Kac formula as 

Zx,.A^,t)= / I?[x(r)] exp - / dr l-^{d^x{r)f - ^^r^{x{r),T) . (2.4) 
Jx{0)=0 \ Jo w 

The average is over all Brownian motion paths starting at at time t and ending at x 

at time t. Thus Z{x, t) can be viewed as the partition function of a continuum directed 

polymer model in 1 + 1 dimensions immersed in a white noise random potential. The 

precise definition of (12. 4 p requires some discussion. We refer to [15] for details. After 

energy renormalization the average partition function is given by 

1 

From the scale invariance of Brownian motion and of white noise one concludes 

that 

Zx,u,d{x, t) = fi ^1,1/2,1 (yU X, p t) (2.6) 

with the notations of (II. 5p . Thus it suffices to study a particular set of parameters and, 
from now on, we set 

A = l, z^ = i, D = l and Zi^i/2,i{x,t) = Z{x,t) . (2.7) 

More important in our context is the observation that the ratio Z{x, t) / {Z{x, t)) for fixed 
t is stationary stochastic process in x. In particular the n-point function is invariant 
after shifting all arguments by the same fixed amount. On the level of the KPZ equation 
this means that under the parabohc shift, h{x,t) +x^/2t, the statistical fiuctuations of 
the interface become stationary in x. 

In the replica approach the n-th moment of Z is written as a matrix element of the 
propagator for the attractive 5-Bose gas. We introduce the Lieb-Liniger Hamiltonian 

^ n ^ n 

Hn = -- - 2 E ^(^^ - • (2-8) 

Then 

(Z(xi, t) . . . Z{xn, t)) = (0, . . . , 0|e-^"*|xi, . . . , x„) , (2.9) 

where \xi, . . . ,Xn) is the quantum state with the particles at positions 

Since the initial condition is a symmetric function, the propagator is needed only in 
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the symmetric subspace. On the level of the 5-Bose gas, the energy renormalization 
corresponds precisely to the omission of the self-energy in (12 .Sp . 
The quantity of interest is the generating function 

n 

G{s-x,t) = (exp (-5^e-^+«(--*))) , (2.10) 

i=i 

setting s = (si, . . . , s„) and x = (xi, . . . , Xn) with xi < . . . < x„. Note that according 
to (ll.4p the term —^x"^ — j^t has been subtracted from h{x, t) already. To compute the 
generating function, one expands the exponential in (I2.10p . writes the multinomials in 
Z(t) as in (12. 9p . and expands the propagator in terms of the exact eigenstates of Hn 
[T3] . As in [14] , for n > 2 this program can be carried out under a natural factorization 
assumption. To distinguish from the exact generating function, we add the superscript 
" whenever the factorization is invoked. 

Our main result is a Fredholm determinant expression for G^. Setting a = (t/2)^/^, 
we obtain 

G\s;x,t) = det{l~QK) , (2.11) 
where K is the Airy kernel (I1.12p and the integral kernel of Q is defined by 

/oo 
du2... dUn (wi|e«^-^^)/2-')^|M2) . . . (M„|e«^"-^i)/2a^)^|w„+i) 
-oo 

X ^{aui - Si, . . . ,aun- Sn) , (2.12) 

with 

pMl _|_ _|_ pUn 

•^("1' • • • , «n) = ' u/' 1 u ■ (2-13) 

1 + e"i + . . . + e^" 

In the long time limit, the shifted height ^{x,t) scales as t^^^, and the joint 
distribution of the ^(x^,t)'s is nontrivial if the distances are scaled as t^/^. Let us 
thus substitute se by aae and xe by 2a^?/^. The right hand side of (I2.10p converges to 
the probability that ^{xj,t) is smaller than aaj for all j, while $ in (I2.12p gives the 
characteristic function of the set {uj > aj,j = 1, . . . ,n}. Hence 

lim Prob(e(2aV, t) < aai, . . . , e(2a'y„, t) < aa„) = F„(a, y) (2.14) 

t—^oo 

and 

F4a, y) = dei{l-K + P.^e^^i-^^^^P^.e^^^-^^)^ . . . P,„e(^"-^i)^i^) . (2.15) 

Here Pa = 1 — Pa is the projector on (— oo,a], in other words {u\Pa\v) = l{u<a}S{u — v). 
The agreement of (12.150 with the standard definition (11.140 will be established in 



Appendix A 



3. Replica Summation 

We start from the generating function G of Eq. (I2.10p . set 

n 

Gi{s;x,t) = (exp(-^e(*/24)-*Z(x^,t))) = G(5;aj,t) , (3.1) 
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where S£ = S£ + x1/2t, and expand the exponential as 

Gi = l + Y.-^l E e-(^^i+-+^^^)(Z(x,„t)...Z(x,,,t)) . (3.2) 

N=l ' fi,...,£iv = l 

By (12.91) the A^-point function of Z can be written as a sum over the eigenstates of the 
(5-Bose gas with particles, 

(Z(7/i, t) . . . Z{yN, t)) = Y. ""''^^ • • • ' ^^1^-) ' (3.3) 



which inserting in (13. 2 p yields 



NtN /2A 



N=l 



^e-*^- 1^.(0,..., 0)p 



^-(s,^+...+s,.,) V'r(a^£i,---,a;^^) 
^ ^.(0,...,0) 



(3.4) 



X > e 

^i,...,^iv=l 

The required matrix elements are computed in [T3l [T6| [T7] . We follow mostly the 
notation in [Hj. Let us pick a number of clusters of particles M, 1 < M < N, and 
the positive integers Ua, a = 1, . . . , M counting the number of particles in each cluster, 
such that 

M 

= iV. (3.5) 

We also introduce M real momenta q^, a = 1, ■ ■ ■ , M, and set g = (gi, . . . , qM), n = 
(rii, . . . , tt-m), y = {Uii ■ ■ ■ 1 Un)- Then the eigenvectors of the iV-particle Hamiltonian 
if AT are indexed by M, q, n and 

see [H], Eq. (3.17). P is the permutation group with elements. The prime means 
that the summation is only over the subset of V consisting of permutations between 
distinct clusters. The coefficients Ap{y) and the sets ^a{p) are defined in [Hj, Section 
III.B, where for the present purpose one only needs to know that 

E'^p(2/) = ^£^0) (3.7) 

and that the sets ^^a(p), a = 1, . . . , M form a partition of the set of integers 1, . . . , A^. 
We now write explicitly the sum over the £a, a = 1, . . . ,N, in (13. 4p . 



E e-K+-+^..)^W(a;,,, . . . = A,{x,,, . . . ,a;,Je 

fl,...,^jv=l P6P 

with £ = {ii, . . . , i^)- The phase 4>{i,p) equals 

M r- n n 

, p) = E - E + E -i 

a=l L j=l j=l 



(3.8) 



1 " 

4 E ~ ^k\ma,j{£,p)ma,k{£,p) 



j,k=i 



(3.9) 
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with 

'ma,M,P) =H-{ce^a{.P)Ac= 3} , j = l,---,'^- (3.10) 

As in the simpler case of the two-point function [H], the sum over p and the 
sum over i are coupled and it does not seem possible to simplify our expression. In 
order to push forward the calculation, we now assume that these two sums factorize. 
While this factorization is not exact, the fact that we recover the n-point distribution 
of the Airy process strongly suggests that the factorization assumption becomes a valid 
approximation in the long time limit. Thus for fixed p we sum the exponential of the 
phase over This sum does not depend on p and (13.71) can be employed to arrive at 

^e-(^^i+-+^^iv)^W(x,„ . . . ,a;,J/^W(0) 



M p n 1 

W exp ^(-Sj + iqaXj)ma,j{t) " 7 XI 1^^' ~ ^k\ma,j{i)ma,k{^) 



e a=l I- j=l 

We use the identity 



j,k=i 



^au+bv+cuv ^cdadt ^au+bv 



which can be proved by expanding both sides as formal series in c, and obtain 



M 



e. a=l 



1] n + kaXj)ma^j 



1 " 

(€) - - X \Xj - Xk\ma,j{i)ma,k{i) 



j,k=l 



M 



-J 12 \^j-^k\ds^dsf. ji; (-Sj+igaa;j)m„,j(£) 

e •''*=^ e-'=^ 



En 

-i E \xj-Xk\dsjds^^ 1 

e 

a=l 



(3.11) 



(3.12) 



(3.13) 



4. ra-point generating function 
4.1. Fredholm determinant 



The generating function Gi with factorization (13. lip will be denoted by G\. The required 
normalization of the eigenfunctions are computed in [161 [13 US] • According to [T^ , Eq. 
(4.5), it holds 

lCHO'---'0)l' = ^'detf^^— ^1-^- . (4.1) 



, 2(n„ + Up) + i(g„ - g/?)/ ^,;3=i,...,M 
In other contexts, determinantal formulas for such normalizations have also been found 
[T8l [To] . The energy of the eigenstate indexed by M, q, n is given by 

^ M M 

EltS = \Y.^^^l-hY.^<-^^)^ (4-2) 



a=l 



OL = l 
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and the properly normalized sum over the eigenstates reads 

M 

\ ( r" \^ 

M 



Ml 

r M=l a=l 



(4.3) 



Then G\ rewrites as 



M=l ■ a=l V n,=l ^ V 2 



:,/3=l,...,Af 



M 



X l[ exp ((tn3/24) - i ^ |x, - Xk\ds^d^ 

a=l j,k=l 

One observes that G\ is equal to a Fredholm determinant 

G| = det(l + R) , 
where R has the kernel 
R{q, m; q' , m') 



(4.4) 



(4.5) 



exp ( {tm^/2A) -\ J2 l^j - Xk\ds^ds^ 

1 ^ i,k=l 



— e 2i 



J2 Q^^jQ-Sj 



27r 



;(m + m') + i(g — q') 



(4.6) 



compare with Eq. (4.9) of [T^ in the two-point case. 

We want to perform the integration over the qa and the summation over the Uq 
inside the Fredholm determinant. For this we first use 



i(m + m') + i(g — q') 



i(m+m')+i((?— g')] 



Then the kernel R can be written as a product of two operators, 

/OO 
dz Ri{q, m; z)R2{z; q', m') , 
-OO 



(4.7) 



(4i 



with 



1 " 

Ri{q,m; z) = t{z>o} exp (^-iqz + {tm^ /2A) - ^ E '^J ~ ^k\ds,ds^ 



and 



i?2(^;g',m') = l{.>o}^e-i'"'V'''^ 



Since det(l + R1R2) = det(l + R2R1), one arrives at 
G\ = det(l + AT) , 

where 



/OO 
dq ^ 7^2(2;; q, m)Ri{q, m; z') . 
°° m=l 



(4.9) 

(4.10) 
(4.11) 

(4.12) 
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More explicitly, it holds 

1 

N{z, z') = l{^,^'>o} exp - I ^ \xj - Xk\dsjds^^ 



j,k=i 

oo 



X 



We now use the identity 

/■oo 



(4.13) 



e 



/oo 
diiAi(ii)e(*/')'^'"^" , (4.14) 
-oo 

in order to perform the summation over m. Defining the function 

gtl _|_ _|_ gWn 

^(ui,...,Un)^- — , (4.15) 

^ ^ 1 + e«i + . . . + e"" ^ ^ 

one has 

N{z, z') ^ -l{.,..>o} exp ( - ^ 1^^- - ^^\9sA>) ^Ai(xi)e^'^(^-^') 

X $ ((t/S)^/^^ - i(^ + ;2') - ^^72 + '^x^q - Sj, j = l,...,n) . (4.16) 

Let us introduce the parameter a — (i/2)^/^. We remove the terms ixjq from $ by 
using 

fiz + a)^exp[ad,]f{z) (4.17) 
and shift uhy u^u + 1^1"^ o?c^ + 2-^/^q;~^(2; + 2;') to obtain 

N{z, z') = -t{z,z'>o} exp ( - ^ \xj - Xk\^,.^,^ [ ^^Ai{u + 

^ ^ j,k=i 

n 

+ 2-^/^a-\z + z')) exp (iq{z - z' 2;^^,.))$ - sj, j^l,...,n). (4.18) 

Integrating over q, 

= 2-V3a-i/2Ai (2-2/3 (6 + a-'/^c)) Ai (2-2/3 (6 - a-'/'c)) , (4.19) 
one arrives at 

/ 1 " /-oo 

Ar(^, ^') = -l{,,,'>o} exp f - - ^ \xj - Xk\ds^d,A / du2-''^/^a~^ 

j,k=i -^-^ 

n n 

X Ai(2-^/\ + a-^z - 2-ia-^ J]x,a3,)Ai(^2-2/3M + a-^z' + 2-^a-i ^ a;,^,.) 
X $ (2-2/3Q,^i _ s^.^ j = 1, . . . , n) . (4.20) 
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Applying (14.1 7p to extract the operator 2^^a^^ X]"=i •^i'^sj from the Airy functions and 
performing the changes of variables u — > and z — az, we obtain 

1 " 

N{z, z) = -aN{az, az') = l{^,^'>o} ^xp - - ^ \xj - Xk\ds^ds^ (4.21) 

j,k=i 

" POO 

— {2a)^^{dz — dz')''^^XjdsA / du Ai(u + z)Ai(u + z')^{au — si, . . . , au — Sn) , 

compare with Eq. (4.19) of [14J. With these rearrangements the generating function Gi 
is written as 

G[ = det{l-N). (4.22) 
We return to the original generating function (under the factorization assumption) 
G^s; X, t) = G{ {s, - Jjx?, ...,Sn- ^,xl; x, t) . (4.23) 
Applying (I4.17P to subtract the parabolic profile, we arrive at 

= det(l - L) (4.24) 

with 



L{z, z') = l{z,z'>o} exp ( - - ^ \xj - Xk\ds^ds^ ^ 5Z ^P'^ ~ (2")"^(5^ - d^') 

j,k=l j=l 

POO 

X ^Xj(9sJ / du Ai{u + z)Ai{u + z')^{au — Si, . . . ,au — Sn) ■ (4.25) 

j=l J-oo 



4-2. Rewriting of the kernel L in terms of the Airy Hamiltonian 

We will now rewrite the kernel L in terms of the Airy Hamiltonian H. In order to 
eliminate the derivatives with respect to the Sj, we first write the function <l> in Fourier 
form as 

n 

/°° -iJ2 rj{au-Sj) „ 

dri...dr„e $(ri, . . . , r„) . (4.26) 

-oo 

Then, setting r = (ri, . . . , r„), the kernel L takes the form 

/oo / \ 

dri ... dvn exp ( i ^ rjSj j L(r)l>(r) , (4.27) 



where 



j,k=l j=l 

\ POO / " \ 

x^Xj-rj] / du exp ( — iau'^^rj \Ai{u + z)Ai{u + z') . (4.28) 



Since the rest of the argument is independent of the particular function $, we will work 
only on L. We use f l4.17p for the derivatives with respect to z and z' and shift u by 
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i{2a)~^ Sj=i ■"^j'^'j integral. Taking into account the fact that Xi < . . . < x„, one 

obtains 

L{r) = l{,,,'>o} cxp (^YjXkrjrk^ exp (^^J^xjrj^ cxp - -^Y^x^j^ 

X J du exp ^ — iau r^^ Ai(M + z)Ai ^ + z' + — 2;jrj^ . (4.29) 

We introduce the operators H, U and D, defined by 

H = -dl+u, U = u, D = d^, (4.30) 

and the function Aiz{u) = Ai{u + z), which is the eigenfunction of the Airy Hamihonian 
H with eigenvalue —z. The Airy kernel K is the projection on the negative eigenstates 
of H, i.e. K\Aiz) = l{z>o}\Aiz). Thus, L{r) can be written as 

L{r) = exp J] Xkrjr-k^ exp Q Xjrf^ ^xp - ^ x^j^ 

X {Aiz'\K exp ( - ^ XjrjD ] exp ( - ia ^ rjU j | Ai^) . (4.31) 

This expression for Lir) will now be rearranged using the following lemma. 

Lemma 1. Let Aj > 0, j = 1, . . . , n + 1 and \Xk G iM, k = 1, . . . ,n. Then the following 
relation holds: 

^XiH^^U^X2H^„,U gA„ffgMnC/eA„+iJf ^ ^^^^ + . . . + A^) 

X exp ( ^(Ai + . . . + Aj)/i| j exp ( 2 HjHk{>^i + . . . + A^) j 

^ j=l ^ ^ l<j<fe<n ^ 

X exp ( - 2 ^(Ai + . . . + Aj)/XjD J exp ( ^ [XjU j 

X exp (^(Ai + . . . + A„+i)/7) . (4.32) 

Proof. The operators C/ and D satisfy the commutation relations 

\E, [/] = -2D , [//, L>] = -1 , [[/, D] = -1 . (4.33) 

Then, the Baker-Campbell-Hausdorff formula yields the commutation relations for the 
exponentials. One has 

Using iteratively these relations, the proof of the lemma is straightforward. □ 
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Lemma [T] with \j = — Xj)/2a'^ and fij = —iarj {xq = Xn+i = 0) allows to 

rewrite the expression f l4.3ip of L{r) as 

L{r) = (Ai^,|i^e("^^/2"')-^e-'""^^e((^i-^2)/2°')^e-'"''2^ 

X . . . e'^^^"-^-^"')/'^"^')^ Q-^'^^'^'^ e^^"/'^'^'^^^ K\Aiz) . (4.37) 

Inserting in f l4.27p . one finds the following expression for the kernel L 

/oo 
du^ ... dUn (Ai,,|ire(-"i/2"')^|Mi)(Mi|e«"^-"^)/2a2)^|M2) x . . . 
-oo 

X (n„_i|e«^"--^")/2a2)^|M„)(M„|e(^"/2a^)^ir|Ai,)<l>(aMi-si^ . (4.38) 

A similarity transformation plus the fact that is a projection, = K, finally gives 
the result announced in Eq. (12. lip . 

5. Conclusions 

Following the strategy in [H] , we extended our results to the spatial n-point generating 
function at some common time t. We also simplified the algebra used for the 2-point 
function. 

The decoupling assumption remains somewhat mysterious. For small M we checked 
that (13. lip is not an identity. Of course, there could be cancellations such that (I4.24p 
survives as an identity for G. We have no idea how to check such a claim directly. On the 
other hand, the long time limit of G" agrees with the conjecture based on universality, 
which is counterintuitive since usually approximations become amplified as t — ?■ oo. 

Our results hold only for sharp wedge initial data. There are two further initial 
conditions for which one would like to solve the KPZ equation. One is the fiat initial 
condition, h{x, 0) = 0. Then 

h{0,t) = log J dxZ{x,t) . (5.1) 

The second one is stationary initial condition. This corresponds to solving ( 12. 2p with 
two-sided Brownian motion as initial condition, i.e. 

Z(±x,0) = e^±(^) , x>0, (5.2) 

with b+ and 6_ two independent standard Brownian motions. Even at the level of the 
one-point function, both cases remain as a challenge. 



Appendix A. Matrix kernel and scalar kernel 

For better readability we use sans-serif for operators on L^(]R) C", while the matrix 
elements, as K^, i,j = 1, . . . ,n, are operators on L^(]R). The operator K of Eq. (11.151) 
can be written as 

K = P(T^K° + T+(K°-1))P, (A.l) 
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where 

Ko-^ijK, P,,- = , (A.2) 

and T", 7+ are lower triangular, resp. strictly upper triangular, according to 

T.^- = l{.>,}c('"-"^)^ , T+. = l{,<,|e(--'"^)^ . (A.3) 
Hence the Fredholm determinant of the operator K can be rewritten as 

det(l - K) = det(l + PT+) det(l - (1 + PT+)-^P(T- + T+)K) . (A.4) 

Recall that 

[H,K] = Q, K' = K, Pl = Pa. (A.5) 

For any integer m > 1, tr(PT+)™ = 0, which implies that det(l + PT+) = 1. Since PT"*" 
is nilpotent, the expansion of (1 + PT"*")"^ terminates as 

(1 + PT+)-^ = 1 - PT+ + . . . + (-1)"-^(PT+)"-^ . (A.6) 

We now consider the operator 

B = (1 + T+)(l + PT+)-^P(T- + T+)(l + T+)-^K'^ . (A.7) 

Since 1 + T+ and K° commute, one has det(l — K) = det(l — B). The inverse of 1 + 7+ 
is given by 

((1 + J+r % = %1 - %„ie("'-"'^+i)^ . (A.8) 
Multiplying on the left by T~ + 7+ leads to 

((T- + T+)(l + J+)-%j = 5y e^"''-"'^)^ . (A.9) 
Finally, setting Pa — l — Pa, one arrives at 

Thus det(l — K) = det(l — Bn), which is the claimed identity. 
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